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SUMMARY 


The vortex-filament method was applied to the simulation of vortex breakdown. The principal vortex 
region was represented by multiple filaments, and an axial velocity component was induced by a spiral wind- 
ing of the filaments. First, an accuracy check was performed for a cylindrical swirling flow with simple 
analytical expressions for the axial and theta velocities. The result suggests that the flow field can be 
simulated to any accuracy by increasing the number of filaments. Second, an axi symmetric- type vortex 
breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial- and 
theta-velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, 
followed axially by a recovery zone and then a second breakdown. When three-dimensional initial data are 
used the second breakdown appears to be of the spiral type in correspondence with experimental observa- 
tions. The present method can easily be used to simulate other types of vortex breakdown or other vortex 
flows with axial velocity. 


1. INTRODUCTION 

The phenomenon of vortex breakdown is not only of fundamental interest, but is of considerable impor- 
tance in aerodynamics as well. For example, vortex breakdown over a wing usually affects the aerodynamic 
performance of that wing. In fact, one of the first experimental studies of vortex breakdown was Peckham 
and Atkinson's investigation (Ref. 1) of flow over a delta wing. At high incidence they observed break- 
down (usually defined as an abrupt change along the vortex axis with a limited region of flow reversal ) of 
the pair of vortices formed at the leading edge of a delta wing. Since then many other experiments of 
vortex breakdown have been carried out. Those studies revealed three basic patterns of breakdown: axi- 

symmetric, spiral, and double-helix. Flow visualizations by Sarpkaya (Ref. 2) show these patterns in 
water. He systematically measured the relationship between Reynolds number, the position of breakdown, 
the pattern of breakdown, and the amount of swirling. Lambourne and Bryer's study of flow over a delta 
wing (Ref. 3; see also Ref. 4) clearly shows two different patterns of vortex breakdown, axi symmetric and 
spiral, occurring at the same time. Recently, Faler and Leibovich proposed a more detailed classification 
scheme for vortex breakdown based on extensive observations (Ref. 5) and obtained measurements of velocity 
profiles using a laser Doppler velocimeter (Ref. 6). Escudier and Zehnder (Ref. 7) have proposed a simple 
criterion for vortex breakdown, using extensive flow-visualization data. See Hall (Ref. 4) and Leibovich 
(Ref. 8) for a review of these studies, as well as theoretical efforts on vortex breakdown. 

On the theoretical side, Benjamin's finite transition theory (Ref. 9) proposes that the breakdown of 
swirling flow corresponds to the transition from supercritical to subcritical flow in a hydraulic jump in 
open-channel flow. On the other hand, Ludwieg (Ref. 10) proposed an instability theory for swirling flows 
in which vortex breakdown is attributed to the growth of instability waves on the vortex core by increas- 
ing helical interference vortices. 

In an attempt to treat vortex breakdown numerically. Hall (Ref. 11) derived a boundary-layer-like 
equation and solved it by a marching method. This scheme assumes axi symmetric, steady flow and can treat 
an arbitrary upstream flow field. Since the equation itself is parabolic, there is some dispute over 
whether this type of equation is valid for sudden changes in flow structure. More recently, Grabowski and 
Berger (Ref. 12) solved the Navier-Stokes equations, using Chorin's artificial compressibility method, to 
simulate axi symmetric vortex breakdown. This scheme also assumes steady flow and appears to be limited to 
low Reynolds numbers. 

It is clearly desirable to attack the vortex-breakdown problem numerically with a three-dimensional 
capability. According to detailed observations, even an axi symmetric- type vortex breakdown includes 
three-dimensional motion. In this paper, we present the result of numerical simulations of vortex break- 
down, using the three-dimensional vortex-filament method developed by Leonard (Ref. 13). In particular, 
we have studied the axi symmetric breakdown. Our present results are for both axi symmetric and three- 
dimensional initial data. It is our view that the primary role of Reynolds number is in shaping the 
internal structure of the vortex core upstream of the breakdown. As shown below, the upstream conditions 
can be duplicated by proper choice of the filament geometry and strength. The process of breakdown occurs 
rapidly and is therefore dominated by nonlinear inviscid phenomena, well described by the dynamics of the 
vortex-filament method. In a similar approach, Del Prete (Ref. 14) modeled vortex breakdown using the 
vortex method. However, she used only parallel longitudinal vortex filaments upstream — with no axial 
velocity component - interacting with a vortex ring to produce a structure suggestive of vortex breakdown. 

Although this report is concerned with the simulation of vortex breakdown, the method used can treat 
any system of vortices with or without axial velocity distribution and, for example, should be capable of 
simulating in detail the merging process of aircraft wake vortices. 


2. NUMERICAL METHOD 

The vortex-filament method was developed by Leonard (Ref. 13) and successfully applied to the wake- 
vortex interaction problem behind large jet aircraft (e.g., B-747) and to the spot-like disturbance in a 
laminar boundary layer. This method seems well suited to the simulation of vortex breakdown, in that com- 
putational elements are required only for the vortici ty-containing fluid in and near the flow region where 
vortex breakdown is occurring. Furthermore, in the vortex filament method rapid changes in the flow 
pattern with three-dimensionality can be represented with minimal error due to numerical diffusion. The 
vortex-filament method is generalized from the two-dimensional vortex-blob method to three dimensions by 
assuming that the vorticity field is represented by a collection of tubes or filaments of vorticity defined 






by the space curves [n(c,t), 1 * 1, 2 L], where £ is a parameter along the curve and L is the 

number of the vortex filaments. The vortex-filament method does not have a singular vorticity distribution 
inside the vortex filament as does the classical point- vortex method. 


The dynamics of these filaments, consistent with the theorems of Helmholtz and Kelvin, is given by 
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The "s" function included in this equation is related to the vorticity distribution inside a vortex 
filament. Leonard (Ref. 13) gave two particular choices of the function "s." In this paper we use the 
so-called "scheme C," which is of the form 


’ 3 [1 + a(o? + a?)/2y 2 ] 3 /^ 

where the choice a = 0.413 for a Gaussian distribution of vorticity within the core yields the correct 
speed of a single vortex ring. 

In vortex-flow simulations, the following two flow representations are possible: 

1. Single filament to represent each physical vortex 

2. Multiple filaments to represent each physical vortex 

The first method was used for test problems which are not described here. These are the interaction 
of four vortex rings and the interaction of two pairs of wake vortices that cross at a right angle with 
some clearance. The results were successful and convincing. In particular, the first case produced 
results very close to those of experiment (Ref. 15). These preliminary results were encouraging for the 
application of this method to vortex breakdown. However, the single-filament representation is appropri- 
ate only when the internal structure is not of major importance. The multiple-filaments representation 
has to be considered when the internal structure, including axial flow, has a significant effect on the 
dynamics. In this paper we use the latter representation to simulate vortex breakdown. If this method 
proves successful and if a detailed procedure is established, many new applications to practical flows 
dominated by vortices may be attempted. 


3. ACCURACY CHECK OF THE VORTEX-FILAMENT METHOD 

We first performed an accuracy check for the three-dimensional vortex-filament method. A flow with a 
simple analytical expression was chosen, and its vorticity was represented by vortex filaments. Then the 
velocity field induced by the filaments was compared with the idealized field. This method has two numeri- 
cal parameters: the core radius and the initial spacing between the vortex filaments. The accuracy is 

expected to improve when the number of the filaments used is increased, which makes the average distance 
between neighboring vortex filaments smaller. It is also interesting to estimate the relation between the 

optimum core radius and the average distance. We did a similar study in two dimensions, with Gaussian 

cores (see Ref. 16). The results suggested the relation o - 6 0 * 875 , where o is the core radius and 6 

is the average distance. This means that the core radius should tend to zero more slowly than the average 

distance does. 

For the test, it is convenient to choose a flow in which the vorticity sharply declines with the dis- 
tance from the axis. The following expressions, often used for unbounded swirling flow, satisfy these 
requi rements: 


V x = a i + a 2 exp(-a 3 r 2 ) (3) 

v e = r„[l - exp(-b 1 r 2 )]/r (4) 

V r = 0 (5) 

where V x is the axial -velocity component, V Q the theta velocity, and V r the radial velocity. The 
values actually used for the parameters are 

a, = 0.95 , a 2 = 1.55 , a 3 = 10 , b 1 = 10.9 , and r 0 = 0.825 

which approximate those obtained in an experiment inside a pipe (private communication, Uchida et al., 
Nagoya Univ.), as shown in Fig. 1. We are not concerned with the wall region where, of course, the 
idealized velocity does not match the experimental velocity because we are simulating an unbounded flow. 
The measurements were obtained by a laser Doppler velocimeter for axisymmetric vortex breakdown at a 
Reynolds number of 2300. Figure la shows the axial velocity component of Eq. (3) which is high near the 
axis and approaches 0.95 (potential flow) as the radius increases. Needless to say, this profile includes 
the axial flow inside the vortex. Figure lb shows the theta velocity, with almost rigid rotation near the 
axis and a maximum value at r = 0.34. Beyond this point the profile becomes proportional to 1/r and 
represents a potential swirling flow for larger values of the radius. Figure 2 shows the theta-vorticity 
distribution, which was calculated from the axial velocity by the relation = -dV x /dr, assuming a 
cylindrical flow. At r = 0.22, the theta vorticity has a maximum value. Figure 3 shows the axial- 
vorticity distribution, which was calculated from the theta velocity by the relation u> x = d(rV e )/dr/r, 
assuming a cylindrical flow. This is a Gaussian distribution. Figure 4 shows the distribution of the 
circulation, which is calculated by the expression 
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This has a familiar pattern of positive curvature near the x-axis approaching a constant value at infinity. 
Figure 5 shows the magnitude of the vorticity vector as a function of r. The deflection angle of the 
vorticity from the x-axis is also shown (dashed line); it increases monotonical ly with radius. Figure 6 
shows a three-dimensional picture of the vorticity vector along the radial direction. It can be seen that 
the vorticity maintains a moderate deflection angle for all radii. 

We now discuss the problem of discretizing the idealized vorticity field with helical vortex filaments. 
The circulation and helicity of each filament as well as their distribution in space must be chosen to 
represent the idealized vorticity field as well as possible. We divided the circular cross section 
radially into a number of zones or rings. Each radial zone in turn contains a number of filaments dis- 
tributed uniformly in the azimuthal direction. In the present simulation we used more filaments per unit 
area near r = 0 where the vorticity is maximum. Using a nonuniform radial distribution of filaments 
allowed us to obtain a more accurate representation of the idealized flow. 

First, we discuss the method used to obtain the filament circulations. For simplicity of analysis, 
we assume infinite, straight vortex filaments with axes in the x-di recti on, distributed continuously along 
the circumference of each ring at the radii r-j . We have the following expression for the axial vorticity 
produced by one of these cylindrical shells of vorticity 


= rx i J J Yri de dx ' (7) 

-oo *0 


where r Xi is the circulation per unit circumference, and y is the distribution function of vorticity 

inside a filament which is calculated from Eq. (2) as follows: 

y(y) = (ds/dy)/4Try 2 = (3aa 2 /4ir)[y 2 + aa 2 ]" 5 ^ 2 
Equation (7) with substitution of Eq. (8) may be integrated analytically to obtain 

2ao 2 r^(r 2 + r 2 + a a?) 

^ 5 = Fxi {[(r-j + r ) 2 

We determine the optimum circulation distribution r Xi for given axial vorticity distribution by minimizing 

the square error with the total circulation fixed. This leads to the following expression for the H 
function 


( 8 ) 
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where C is the total circulation. Differentiating H with respect to the r Xi and X yields the follow- 
ing simultaneous equations for r Xi : 
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Solving Eqs. (11), we obtain the optimum circulation distribution r Xl - for each ring. If N-j is the 
number of filaments placed in the ith ring then the circulation of each of these filaments is given by 

T j = 2Trr i r x ./Ni (12) 


The next step is to determine the angle each helical filament makes with the x-axis, the deviation 
angle. For this purpose we assume that the theta component of vorticity is distributed continuously and 
uniformly along the x-direction with r 6l - circulation per unit axial length for the ith ring. Again, 

using Eq. (8) we find that the ith ring of distributed theta vorticity contributes to the theta vorticity 
component as follows 


^e-j 


1 [(rf 


4zr 2 aa| 

r 2 ) 2 + 2aaf(r? + r 2 ) + a 2 a]V/ 2 


(13) 


The same minimization procedure used above was used to determine the optimum except that total vor- 
ticity constraint is not imposed in this case. Finally the deviation angle from the x-axis of each fila- 
ment in ring i was determined by 


tan oj = r 9i /r Xi 


(14) 
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By using the circulations and deviation angles for each ring, we calculated the induced velocity at a 
cross section using the following expressions: 


L 0 00 

= -(l/4w) 2 rj J ajRj[r sin(ajx' + 0j) - Rj ]/F dx' 
j=: -oo 


Vy * E (Vy sin §j - V z cos 0j) 


L 

V z * (Vy cos §j + V z sin e j ) 




(15) 


where 


F = [x ' 2 - 2yRj cos(ajx') - 2zR^ sin(ajx') + Rj + aoj + r 2 ] 3 / 2 


- ( 1/4tt ) r j J [z sin 0j - y cos §j - Rj sin(ajx') + ajRjX 1 cos (a jx ' )]/F dx' 

- ( l/4ir)rj r [-y sin 0j - z cos 0j + Rj cos(ajx') + ajRjX* si n (a jx ' ) ]/F dx' 


Pj = 2 tt/Xj , e j = e j “ 6 

(Rj,0j) and (r,0) are the coordinates of the filament and the point where the induced velocity is required, 
respectively, and pj is the filament pitch. These expressions are obtained by transforming the original 
expressions of the induced velocities, written in the Cartesian coordinates, using the relations repre- 
senting the spiral shape of the filaments. 


As our first test of our discretized vorticity field we calculated the distribution of the induced, 
unperturbed axial and theta velocity components of the upstream flow, which for the test extended over the 
entire x-axis. The actual integration was performed numerically from -10 to 10 using Simpson's rule with 
a spacing of 0.04 for x. Four rings were used radially with radii at 0.06, 0.18, 0.30, and 0.42 with 
6, 6, 8, and 8 filaments in the respective rings. The core radius o of each filament was 0.16. The 
results along the z-axis are depicted in Fig. 7 and show very good agreement. To check for uniformity in 
the azimuthal direction we obtained contour plots of the induced velocity field in a plane normal to the 
x-axis. These are shown in Fig. 8 along with corresponding contours of the idealized field. Generally we 
obtained a very good reDresentation of the idealized field with a tendency to slightly underestimate 
idealized values. 


4. SIMULATION OF VORTEX BREAKDOWN 

For conditions upstream of vortex breakdown we assumed the idealized flow from above and therefore we 

used the filament parameters determined above. There might be some dispute over the difference between the 

experimental flow inside the pipe and the unbounded flow assumed for the simulation. As mentioned above, 
stable vortex breakdowns are produced inside the pipe with a divergent wall or without. Although some 
differences might exist between the two configurations, we assume the essential ingredients remain for 
vortex breakdown in unbounded flow. Experimentally, vortex breakdown is easily observed in a pipe at mod- 
erate Reynolds number which makes quantitative measurements less difficult. 

In its present form the vortex filament method does not explicitly include the effect of viscous dif- 
fusion. Therefore, we cannot say anything about the effect of Reynolds number on the breakdown process 
itself. We assume this process to be dominated by nonlinear inviscid effects. However, the upstream pro- 
files within the vortex depend on Reynolds number as well as the vane angles which control the amount of 
swirling. However, by shaping the vortex filaments upstream of breakdown we can reproduce any desired 
velocity (and hence vorticity) profiles. In this sense we can follow changes in Reynolds number, and in 
the present simulation we attempt to simulate a vortex breakdown occurring at an experimental Re = 2300. 

For numerical purposes, we divided the flow domain into four axial subregions (Fig. 9). For simplicity 

just one filament is shown. Regions I, II, and IV have predetermined filament geometries. Only in 

region III is the flow simulated fully. Regions I, II, and IV are constructed to provide high-quality 
boundary data for region III. Regions I and IV consist of straight, semi -infinite filaments that induce 
a swirl component on the filaments in region III. Here, analytical expressions were used for the induced 
velocity. For example, if one of these semi-infinite filaments, parallel to the x-axis, extends from 
(x 0 , y 0 , z 0 ) to minus infinity, the analytical induced velocity is as follows: 


v y * -(rMn)U - z 0 )F K (16) 

v z = (r/4ir)(y - y„)F J 
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where 


F = { (x 0 - x)/[(x - x 0 ) 2 + (y - y 0 ) 2 + (z - z Q ) 2 + A] 1 / 2 + l}/[(y - y 0 ) 2 + (z - z 0 ) 2 + A] 

A = ot(a 2 + a j )/2 

Region II consists of spiral or helical filaments which rotate and translate with fixed, predetermined 
velocities. The filaments in this region also induce velocity on the filaments in region III. Filament 
geometries and strengths in region II were determined as in the previous section. Thus the filaments in 
region III move according to an induced velocity from the filaments of region III as well as all other 
regions, using the dynamic equation, the Biot-Savart Law. 

A Runge-Kutta fourth-order scheme was used for time-integration. Special attention was given to the 
filaments of region II. Because the Runge-Kutta scheme requires four different substeps to obtain the 
final solution at each time-step, we must use three different geometries other than the "exact geometry" 
to simulate precisely the Runge-Kutta errors that are being made in the upstream part of region III. Both 
the radius and phase angle were precalculated for each substep and stored in the program. 

At each time-step node points were redistributed and, if necessary, added to each filament to maintain 
a nearly constant distance between points. This was accomplished by using a parametric spline fit to 
obtain the space curve defining a given filament. 

At present we have insufficient knowledge about the initial disturbance necessary to produce vortex 
breakdown. Each type of vortex breakdown might need a different initial perturbation. At least it might 
be natural to consider that the three-dimensional type breakdown will need a three-dimensional disturbance. 
Preliminary axi symmetric calculations showed that a small amplitude disturbance on a vortex with no axial 
flow resulted only in weak traveling waves, but a large amplitude disturbance causes breakdown similar in 
appearance to wave-breaking. Experimentally, vortex breakdown is sometimes very sensitive to a disturbance 
in that a change in breakdown pattern occurs even if the same upstream conditions are maintained. In this 
report we present results for two initial disturbances. Case I had an axi symmetric disturbance and Case II 
had a three-dimensional disturbance. 

The top function (t = 0) in Fig. 10a shows the initial geometry of the vortex filaments used for Case I 
and represents a large amplitude disturbance. The region of full simulation (region III) extends over the 
range 0<x<5. In 0<x<3, the helical geometry of the upstream flow was used, but in 3 < x < 4, the 
helices gradually reverse pitch and the pitch was elongated by a factor of 1.4 in this decelerating region. 
In 4 < x < 5, straight filaments were used. 

Figure 10 also shows the time variation in the vortex filaments. The time-step was 0.2. A uniform 
axial velocity of 0.4 (potential flow) was added to that induced by the vortex filaments. This value was 
chosen after the results of preliminary runs were available, giving us an estimate of the speed of the 
disturbance. If this flow component is not chosen with some accuracy the disturbance will rapidly move 
toward one of the computational boundaries. However, this value does not play an essential part in flow 
pattern, because this is just the moving frame through which we see the phenomenon. Alternatively, one 
could continually adjust the free-stream component during the simulation to keep the disturbance centered. 
During the first few time-steps, a remarkable vortex swelling can be seen in the downstream region where 
the initial disturbance had a flow reversal (t = 1). This swelling moves downstream (t = 2). At later 
times these vortex filaments are elongated in the axial direction and in particular the filaments in the 
third ring appear almost parallel to the x-axis (t = 2 and t = 3). Just downstream of this region a 
narrow region with small pitch can be observed (t = 3). A central region which appears devoid of filaments 
due to the elongation and straightening of the filaments can also be seen at t = 6 in the middle of the 
picture. This region moves upstream and is a manifestation of the axi symmetric vortex breakdown. On the 
other hand, at t = 4 the filaments in the second and third rings show disorder just downstream of the 
center of the picture, followed by a contraction toward the x-axis near the downstream boundary. In con- 
trast with this, the filaments in the outermost ring show slight expansion near the downstream boundary 
with some contraction upstream. At later times the outermost filaments establish a wavy pattern downstream 
of the first expansion. At t = 6 most filaments have become very disordered near the downstream boundary. 
Just upstream of this expansion the filaments of the second and third rings have become "untied" and elon- 
gated in the axial direction. Another interesting feature is that the filaments in the innermost ring 
contract to r = 0 at the station where the outer filaments show the expansion at t = 2. This region of 
contraction has a long axial extent by t = 3 and then moves downstream with time. 

In Case II we added a three-dimensional perturbation to the initial disturbance as shown in Fig. 10b 
(t = 0). The perturbation was caused by adding a radial displacement to each filament which was a sinu- 
soidal function of theta (wavenumber = 1) and was nonzero and sinusoidal in x for 3 < x < 4. It can be 
seen that differences in the filament configurations between Case I and Case II persist in time but remain 
confined to the downstream end of the computational domain. 

Figure 11 shows the contours of axial velocity. Solid lines denote positive values and dashed lines 
negative values. Figure 11a shows Case I. Initially (t = 0), a downstream region of low axial velocity is 
shown as expected from the initial geometry of the filaments. By t = 1 this low-velocity region has gone 
through the downstream boundary and cannot be seen. Instead a high-velocity region appears. Between t = 1 
and t = 2 this high-velocity region enlarges and moves downstream. At later times we see the development 
and enlargement of a low-velocity region upstream of the high-velocity region. Meanwhile, the peak velocity 
in the high-speed region is decreasing. At t - 7 and t = 8 we see another interesting feature. Near 
the downstream boundary another low-velocity region appears. Apparently, a second breakdown is occurring 
near the downstream boundary. Experimentally it is very often observed that axisymmetric vortex breakdown 
is followed by a second breakdown of the spiral type. In the experiment used to determine the upstream 
conditions a clear spiral-type vortex breakdown was observed. However, the axisymmetric initial distur- 
bance for Case I only leads to axisymmetric flow. 
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Next, we discuss Case II (Fig. lib), having a three-dimensional initial disturbance. Generally, 

Case II shows a pattern similar to that of Case I. However, we again see that three-dimensionality 
remains at later times, but only near the downstream boundary. As seen at t = 5 and later, the contours 
in the first breakdown are very similar to those in Case I and appear axi symmetric. But the second break- 
down is distorted, not axi symmetric, and appears to be rotating around the x-axis (as seen in the down- 
stream region at t = 7 and 8). This is suggestive of a vortex breakdown of the spiral type. 

Figure 12 shows the theta velocity contours. For Case I at t = 0, the contours show no disturbance 
in the theta velocity because the filaments are distorted initially only in the theta direction. At t = 1, 
swelling of contour lines can be seen in the region of the initial perturbation. At t = 2 an interesting 
phenomenon appears downstream. Another swelling after the first swelling can be seen near the downstream 
boundary. At this time the swelling of the filaments near the x-axis seen at t = 1 has disappeared. 

The contours at t = 3 are similar to those at t = 2. Between t = 2 and t = 4 the axial position of 
the maximum swelling of the level 3 contour clearly moves upstream. After t = 6 the first swelling 
spreads to the contours near the x-axis. The first and the second swellings correspond to the two break- 
downs observed above. 

For Case II (Fig. 12b) we see that at t = 0 the contour lines have a slightly crooked pattern for 
3 < x < 4 where the initial three-dimensional disturbance was imposed. At t = 1 we can see the clear 
three-dimensional wavy pattern inside the region of swelling. At t = 2 the three-dimensionality moves 
toward the downstream boundary and seems to disappear until t = 7. At t = 7 and t = 8 we clearly 
observe three-dimensionality near the downstream boundary and a comparison of the contours at these two 
times indicates a rotating disturbance. This is in contrast with the first breakdown whose contours 
appear to be nearly axi symmetric and agree well with those of Case I. We conclude that the complete pat- 
tern of an axi symmetric vortex breakdown followed by a breakdown of the spiral type was simulated by allow- 
ing for three-dimensionality. 


5. CONCLUDING REMARKS 

The vortex-filament method was applied to the numerical simulation of vortex breakdown, which is 
intrinsic to the longitudinal vortex. Generally, there are three types of vortex breakdown: axi symmetric, 

spiral, and double helix. In this paper we simulated the axisymmetric type of vortex breakdown. Generally, 
it is difficult to obtain information experimentally due to the sensitivity of breakdown to any perturba- 
tion. Calculated results should be helpful to understand this complicated phenomenon. The principal 
vortex region was represented by the multiple filaments. One filament is not enough to represent this 
complicated flow field. In particular, the present method can include axial flow which is often observed 
in physical vortices. 

Vortex breakdown was simulated under the assumption of nonlinear, inviscid dynamics of vorticity using 
the three-dimensional vortex-filament method. Viscous effects are assumed to be important only in forming 
upstream velocity profiles. Upstream profiles were simulated accurately with helical vortex filaments by 
choosing pitch and circulation to minimize the error. First, an accuracy check was performed for a cylindi- 
cal flow with simple analytical functions for the axial and theta velocities. This profile was the idealiza- 
tion for the experimental data used for the vortex-breakdown simulation. Comparison of the simulated 
results with the idealized flow suggests that the flow field can be simulated to any accuracy by increasing 
the number of the filaments. 

Two cases. Case I and Case II, were calculated to simulate axi symmetric -type breakdown. In order to 
create a strong initial disturbance, the filaments in the downstream axial region were provided with a 
change in pitch. We used two kinds of initial disturbance to see their effect on vortex breakdown. Case I 
had an axisymmetric disturbance and Case II had a three-dimensional disturbance. Filament configurations 
and contour plots of velocity showed, for both cases, two simultaneous vortex breakdowns. One vortex break- 
down was similar in both cases and was axisymmetric. The other vortex breakdown was axisymmetric for axi- 
symmetric disturbances and three dimensional for three-dimensional disturbances. The latter case agrees 
with the experimental observations where an axisymmetric vortex breakdown is followed by a spiral-type 
breakdown. 

Other types of breakdown (spiral and double-helix) should be possible to simulate by this method. 

The present method can also be applied to other vortex flows with axial velocity profile. 
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Fig. 1. Velocity profiles in the upstream region, (a) Axial velocity; (b) Theta velocity. 




Fig. 2. Upstream theta vorticity. Fig. 3. Upstream axial vorticity. 
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Fig. 4. Circulation. 


Fig. 5. Norm and angle of the vorticity vector of the 
idealization as a function of radius. 
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AXIAL VELOCITY 




Fig. 7. Comparison between the calculated velocity distribution along the z-axis and the idealization, 
(a) Axial velocity; (b) Theta velocity. 
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Fig. 8. Upstream velocity contours in a cross section, (a) Axial velocity; (b) theta velocity. 
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Fig. 9. Computational regions for filament to calculate induced velocity. 
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Fig. 10. Evolution of vortex filaments, (a) Case I; (b) Case II. 
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Fig. 11. 


Axial velocity contours. 


(a) Case I; (b) Case II. 
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Fig. 12. Theta velocity contours, (a) Case I; (b) Case II. 
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